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I. INTRODUCTION 

The manipulation and control of magnetic domain 
walls (DWs) in ferromagnetic nanowires has recently be- 
come a subject of intense experimental and theoretical 
research. The rapidly growing interest in the physics of 
the DW motion can be mainly explained by a promising 
possibility of using DWs as the basis for next-generation 
memory and logic devicesi"— . However, in order to realize 
such devices in practice it is essential to be able to posi- 
tion individual DWs precisely along magnetic nanowires. 
Generally, this can be achieved by either applying ex- 
ternal magnetic field to the nanowire, or by generating 
pulses of spin-polarized electric current. The current 
study is concerned with the former approach. 

Even though the physics of magnetic DW motion under 
the influence of external magnetic fields has been studied 
for more than half a century-"- , current understanding of 
the problem is far from complete and many new phenom- 
ena have been discovered only recentlyiSr— . In particu- 
lar, a new regime has been reporte d^'^'^^ in which rigid 
profile DWs travel along a thin, cylindrically symmetric 
nanowire with their magnetization orientation precessing 
around the propagation axis. In this paper we address 
the stability of the propagation of such precessing DWs 
with respect to perturbations of the initial magnetization 
profile, some anisotropy properties of the nanowire, and 
applied magnetic field. 

Let m(a;) = (cos6'(a;), sin 6'(a;) cos (/>(a;), sin6'(a;) sin 0(x)) 
denote the magnetization along a one-dimensional wire. 
With easy magnetization axis along x and hard axis 
along y, the micromagnetic energy is given hy^ 



E{m) = 




= i y" (^Ae'^ + sin^ e{A(f>'^ + K1+K2 cos2 (j))^ dx (1) 

where A is the exchange constant and Ki, K2 the 
anistropy constants. Here and in what follows, integrals 
are taken between —00 and 00 (for the sake of brevity, 
limits of integration will be omitted) . 

We consider here the case of uniaxial anisotropy, K2 = 



0. Minimizers of E subject to the boundary conditions 
lim m{x) — ±x, (2) 

X— >±oo 

describe optimal profiles for a domain wall separating 
two magnetic domains with opposite orientation. The 
optimal profiles satisfy the Euler-Lagrange equation 

m X H = 0, (3) 

where 

SE 

H = — - — = Am" + Ki(m ■ x)x = — Com + ein -I- e2p. 
dm 

(4) 

Here m, n = dm/ 89 and p = m x n form an orthonormal 
frame, and the components of H in this frame are given 

by 

eo = A9'^ + sin^ diKi + A(f>''^) 

ei = AO" - i sin2e'(i^i + A(t>'^), 

62 = AsinOcI)" + 2Acosee'(t)'. (5) 

In terms of these components, the energy Eq. ([TJ (with 
K2 = 0) is given by 

^(m) = eodx, (6) 

and the Euler-Lagrange equation becomes ei = 62 = 0. 

While the energy E is invariant under translations 
along and rotations about the x-axis, the optimal pro- 
files cannot be so invariant (because of the boundary 
conditions). Instead, the optimal profiles form a two- 
parameter family obtained by applying translations, de- 
noted T(s), and rotations, denoted i?(cr), to a given op- 
timal profile m*. We denote the family by T{s)R{a)m^,. 
In polar coordinates, T{s)R(a)m^, is given by </'(x) = a 
(the optimal profile lies in a fixed half-plane), and 0{x) — 
9f.{{x — s)/do), where do = \J AjKi and 

6',(0 = 2tan-i(e"«). (7) 
It is clear that 6'*(C) satisfies 

B',^- sin 61, , sin (0 = sech ^ (8) 
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The dynamics of the magnetization in the presence of 
an apphed magnetic field is described by the Landau- 
Lifschitz-Gilbert equationi^, which for convenience we 
write in the equivalent Landau-Lifschitz (LL) form, 

m = m X (H + Ha) - am X (m X (H + Ha)) . (9) 

Here a > is the damping parameter, and we take the 
applied field to lie along x, 

Ha = i/i(t)x. (10) 

In polar coordinates, the LL equation is given by 

6 — aei — 62 — OiHi sin0, (11) 
sin 6*0 = ei + ae2 — Hi sin 6*. (12) 

The precessing solution is a time-dependent translation 
and rotation of an optimal profile, which we write as 
T{xQ{t))R{(t>Q{t))m^. The centre a;o(t) and orientation 
00 {t) of the domain wall for the precessing solution evolve 
according to 



- ado Hi, 4>o 



-Hi. 



(13) 



It was shown-'^'^'-'^'* that T(a;o)i?(0o)m* satisfies the LL 
equation. 

It is important to note that the precessing so- 
lution is fundamentally different from the so-called 
Walker solution®. Indeed, the latter is defined 
only for K2 > (the fully anisotropic case) and 
time-independent Hi less than the breakdown field 
Hw — aK2/2. The Walker solution is given by m(x, t) — 
( cos 6w{x, t), sin 6w{x, t) cos sin 9w{x, t) sin 
with 



ew{x,t) ^ e,{-i-\x -Vwt)) , 

sin20vv — Hi/Hw , 



and 



Vw = lia 



')doHi , 



Ki 



Ki + K2 cos^ 4> 



w 



(14) 
(15) 



(16) 
(17) 



Equations (fT4|) - ([T7)) describe a DW traveling with a 
constant velocity Vw whose magnitude cannot exceed 
7(a -I- a~'^)doHw', note that Vw does not depend lin- 
early on the applied field Hi. In contrast, the velocity ig 
of the precessing solution is proportional to Hi , and can 
be arbitrarily large. Also, while for the Walker solution 
the plane of the DW remains fixed, for the precessing 
solution it rotates about the nanowire at a rate propor- 
tional to Hi. Finally, for the Walker solution, the DW 
profile contracts (7 < 1) in response to the applied field, 
whereas for the precessing solution the DW profile prop- 
agates without distortion. 

In this paper we consider the stability of the precess- 
ing solution. We establish linear stability with respect 



to perturbations of the initial optimal profile (Sec. Hit , 
small hard-axis anisotropy (Sec. IIIip . and small trans- 
verse applied magnetic field (Sec. IIV[) : specifically, we 
show, to leading order in the perturbation parameter, 
that up to translation and rotation, the perturbed solu- 
tion converges to the precessing solution (in the case of 
perturbed initial conditions) or stays close to it for all 
times (for small hard-axis anisotropy and small trans- 
verse magnetic field). The argument is based on consid- 
erations of energy, and depends on the fact that for all 
t, the precessing solution belongs to the family of global 
minimizers. The analytic argument establishes only lin- 
ear stability. Nonlinear stability is verified numerically 
for all three cases in Sec. [V] For convenience we choose 
units so that A ^ Ki — 1. 



II. PERTURBED INITIAL PROFILE 

Let m.^{x,t) denote the solution of the LL equation 
with initial condition a perturbation of an opti- 

mal profile. Let T{x^{t))R{(t)^{t))m^, denote the optimal 
profile which, at time t, is closest to m.^\ that is, the 
quantity 

||me-T(s)i?(cr)m*|p = j {iaf_{x,t)-R{(T)ia^{x-s)Y dx 

(18) 

b^{t). Then the 



is minimized for s — X({t) and a 
following conditions must hold: 



T(x,(t))i?(0,(i)) 



9m* 

dx 



dx — 0, 



/ 



m, • (x X T{xe{t))R{(f>e{t))m^) dx = 0. 



(19) 



It is clear that x^{t) = XQ{t) + 0(e) and 0e(i) = 
4)o{t) + 0{e), but we shall not explicitly calculate the 
0(e) corrections produced by the perturbation. Rather, 
our approach is to show that to leading order O(e^), 
||me — r(a;£)i?((/)e)m*|p decays to zero with t. This will 
imply that the precessing solution is linearly stable under 
perturbations of initial conditions up to translations and 
rotations. 

Let 9g{x,t) and (t>eix,t) denote the spherical coordi- 
nates of mc(a;,t). We expand these in an asymptotic 
series, 

9,{x,t) = e^x - x,{t)) + 66*1 (x - Xe{t),t) + ■■■ , 
(l),{x,t) ^ <j),{t) + e(j)i{x-x,{t),t) + --- (20) 

where the correction terms 6'i(^,t), 0i(^,t), etc are ex- 
pressed in a reference frame moving with the domain wall 
Then to leading order O(e^), 



me - r(a;e)i?(0c)m* 



9? +sin2 6**02) 



h \0i) + e^ (sin 6**01 |sin6'*0i} , (21) 
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where for later convenience we have introduced Dirac no- 
tation, expressing the integral in Eq. (PT|) in terms of 
inner products. It is straightforward to show that the 
conditions Eq. ([T^ imply (using 9'^ = — sinf?*)) that 

(sin 61, \ei) = (sin 6*, |sin6',0i) = 0, (22) 

which expresses the fact that the perturbations described 
by 6i and 0i are orthogonal to infinitesimal translations 
(described by sin0,) along and rotations about x. 

Since the difference between ra^ and T(a;e)i?(0e)m, 
is 0(e), the difference in their energies is O(e^) (as 
T(a;e)i?((/>c)mH. satisfies the Euler-Lagrange equation 
Eq. ([3])), and is given to leading order by the second 
variation of E about m,, 

A£;, = E{m.,) - E{T{x,)R{(t),)xn.^) = 



E{in^) ~ E{m.^) = 




where /o = + cos26i*6i? + sin^ e^^\^ . 

Using the relations Eq. (jSj and performing some integra- 
tions by parts, we can write 

j hd£.^ {Oi \'H\ei) + (sin e,(t)i I U |sin e,<j){) , (24) 

where % is the Schrodinger operator —<P /d^^^V{S) with 
potential given by 

V{i) = 1-2 sech^ (25) 

V{^) is a particular case of the Poschl- Teller potential, 
for which the spectrum oiH is known^^. H has two eigen- 
states, namely sin0,(^) — sechf with eigenvalue Aq = 0, 
and cos6',(^) — tanh^ with eigenvalue Ai = 1, and its 
continuous spectrum is bounded below by A = 1. This 
is consistent with the fact that the optimal profiles are 
global minimizers of E (subject to the boundary condi- 
tions Eq. ([5])), which implies that the second variation 
of E about is positive for variations transverse to 
translations and rotations of m,. It follows that, for any 
(smooth) square- integrable function /(^) orthogonal to 
sin 9^ , we have that 

{f\W+^\f)>{f\W\f) (26) 

for j > (we will make use of this for j — Q and j = 1). In 
particular, since 9i and sin0,0i are orthogonal to sin0* 
(cf Eq. ([22])), it follows that 

> (01 (27) 
(sin6'*(/)i| ?^ |sin6'*(/)i) > (sin0*0i \s\ii9^4)i) . (28) 

Therefore, from the preceding Eqs. (f27 |) -([28 |) and 
Eqs. ^ and ([23])" ([241), we get, to leading order 0{e^), 
that 

||m, -r(x,)i?(0,)m*|p < 2A£;e. (29) 



Below we show that, to leading order O(e^), for small 
enough Hi (it turns out that \Hi\ < 1/2 is sufficient), we 
have the inequality 

j^AE, < -^AE, (30) 

for some 7 > 0. Taking Eq. ([50)1 as given, it follows from 
the Gronwall inequality that 

AE, < \ce^e-^' (31) 

for some C > (which depends only on the form of the 
initial perturbation). From Eq. (|29|) . it follows that 

\\m, - r(x,)i?(0e)m* IP < Ce^e"'^*. (32) 

The result Eq. ([32]) shows that, to O(e^), nie converges 
to an optimal profile with respect to the L^-norm. In fact, 
with a small extension of the argument, we can also show 
that, to O(e^), converges to an optimal profile uni- 
formly (that is, with respect to the L°°-norm). Indeed, 
making use of the preceding estimates, one can obtain 
a bound on ||m^ — T{Xf^)R{(j)^)TLn'^,\\, the L^-norm of the 
difference in the spatial derivatives of the perturbed so- 
lution and the optimal profile. To O(e^), 

||m^-T(a;e)i?(</>.)m'Jp 
= {{9[ \9[) + (sm0^(j)[ |sin6»*0;) + (sin 6*, 6*1 |sin6l,6ii)) 
< (3((6'i|?^ |6li) + {sin0,(f>i\H \sm9,(f>i)) 

< Qe^AE,. (33) 

Arguing as in Eqs. ((29)) "([32 |) . we may conclude that 
||m^ — T{xe)R{(j)^)va'^\\ decays exponentially with t. 
Thus, nig converges to an optimal profile with respect 
to the Sobolev Hi-norm (where = + ||/'|1^). 

It is a standard result that this implies that the conver- 
gence is also uniform (again, to 0(e^)). 

It remains to establish Eq. pop . From Eq. we have 
that for any solution m(a;, t) of the LL equation, 

— E{m) = — / H • riida; 
at J 

— J (m X H) • Ha dx— 
- a J [m X H)^ + (m X H) • (m X Ha) dx 

= -'^ J (e? + el + Hi sinSei) dx, (34) 

where ei and 62 are given by Eq. ([S]), and we have used 
the fact that the term (m x H) • Ha vanishes on inte- 
gration. Substituting the perturbed solution m,; into 
Eq. (1311) and noting that the E{T{x^)R{(j)^)m^) = E{m^) 
does not vary in time, we obtain after some straightfor- 
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ward manipulation that 



dt 

-ae^{{ei\H^ \ei) + (sin 61,011^2 |sin6l*<?!)i) + i/iF ) 

(35) 

to leading O(e^), where 

F = I (cos 0Jo + cos 9^ sin^ 9^91) dS,. (36) 



For the first two terms on the rhs of Eq. (1551) . we have, 
from Eq. ((26]) and Eqs. ([231) -(EH), that 



{0i\'H'' \9i) + {sm9,(j)i\H^ |sin6l*^i) 

> {9i\n\9i) + (sin6',</>i|H|sin6i,(/)i) 

- ^AE,. (37) 

The term HiF in Eq. ([35]) is not necessarily positive, 
as Hi can have arbitrary sign. But for sufficiently 
small I Hi I , it is smaller in magnitude than the preceding 
two terms. Indeed, we have, again using Eq. ([26]) and 
Eqs. ([23])-([ll]), that 

|i^| < J (l/ol + Oi^) d^ < ^AE, + {9i \9i) 



< -^AE, 



h\H\9i) <-AE,. (38) 



Substituting Eqs. ([37]) and ([3H]) into Eq. ([HS]), we get that 
d 



dt 



AE, < -2a(l - 2\Hi\)AE, 



(39) 



from which the required estimate ([30]) follows for \Hi \ < 
1/2. 

It is to be expected that the stability of the pre- 
cessing solution depends on the applied field not being 
too large. Indeed, it is easily shown that, for Hi > 1 
(resp. Hi < —1), the static, uniform solution m = — x 
(resp. m = +x) becomes linearly unstable. As the 
precessing solution is nearly uniform away from the do- 
main wall, one would expect it to be similarly unstable 
for |iJi| > 1. The numerical results of Sec. IV Al bear this 
out. Finally, we remark that the stability criterion ob- 
tained here, namely \Hi \ < 1/2, is certainly not optimal. 



III. SMALL HARD-AXIS ANISOTROPY 

Next we suppose the hard-axis anisotropy is small but 
nonvanishing, taking K2 — e > 0. Let m.^{x,t) denote 
the solution of the LL equation with initial condition 
m^{x,0) — m*(x). As above, let T{xe{t))R{(j)f{t))mf 
denote the translated and rotated optimal profile closest 
to nic at time t. Adapting the argument of the preceding 
section, we show below that, to leading order O(e^), 

||m, - T{xe)R{(t),)ni^\\^ < Cse^ for ah t > (40) 



for some constant C2 > 0. In contrast to the preceding 
result Eq. (|32|) for perturbed initial conditions, here we 
do not expect nic to converge to r(xe)i?(0£)m*. Indeed, 
while an explicit analytic solution of the LL equation is 
not available for small K2 (the Walker solution is valid 
only for K2 > 2\Hi\/a), it is easily verified that there are 
no exact solutions of the form T {x^{t))R{(l)^{t))m^, . The 
result Eq. (j40p demonstrates that, through linear order in 
e, the solution for K2 = e remains close to the precessing 
solution, up to translation and rotation. 

To proceed, let AE^ denote, as above, the difference in 
the uniaxial micromagnetic energy, i.e. the energy given 
by Eq. ([1]) with K2 — 0, between nig and T{xe)R{(j)e)m^:. 
Then, as in Eq. (|29l) . we have that 



|m, - T(a;,)i?(0,)m*|p < 2AE, 



(41) 



As E{T{xe)R{<j)e)^*) ~ E{m^,) is constant in time, we 
have that 



(42) 



The hard-axis anisotropy affects the rate of change of 
the uniaxial energy through additional terms in rh. In- 
deed, for any solution m(x, t) of the LL equation, we have 
that 



dt''^'^^ - dt 



E{m) + G{m), 



(43) 



where d/ dt\K^=oE{\n) denotes the rate of change when 
= 0, as given by Eq. ([M]). and 

G{m) — I (m • y)(m X H(m)) • y dx 

+ ea J {m X H(m)) • (m x y)(m • y) dx. (44) 

Taking m = nig, we recall from the preceding section 
(c.f. Eq. ([20])) that, for \Hi \ < 1/2, 



dt 



E{m,) < -jAE^ 



(45) 



K2=0 



for some 7 > 0. Below we show that there exists con- 
stants Ci , 7i with 7i < 7 such that 



|G(m,)| <7iAS, + Cie2 



(46) 



Taking Eq. (I46p as given and substituting it along with 
Eq. ^ into Eqs. ([12])-([13]), we get that 



-AE, < -{j-^i)AE, + Cie\ 
From Gronwall's equality it follows that 

AE, < -^e\ 
7-71 



(47) 



(48) 
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which together with Eq. (|4ip yields the required result 
Eq. (HOD. 

It remains to show Eq. (|46p. Substituting the asymp- 
totic expansion Eq. ((^0)) , we obtain after straightforward 
calculations that, to leading order O(e^), 

G(me) = -e^cos^ 

X J {sin* 9^(j)[+ A/ 3a sin^ 9 J[) dt (49) 

This can be estimated using the elementary inequality 



2\ab\<f3a' + j, 



(50) 



which holds for any /3 > 0. Indeed, recalling Eqs. ([S]), 
(|23| . (|27)) . and using integration by parts where neces- 
sary, we have that 

J sin* e,<j)[d^ <^ J sin^ e,(P[^ + ^ J sin^e^de 

J sin^ e,e[ d^ < ^ / + 2^ / ^'''^ ^* 



15/3' 



(51) 



From Eqs. (|49)) - ([5T|) . it is clear that /3, 71 and Ci can be 
chosen so that Eq. pS)) is satisfied. 



IV. SMALL TRANSVERSE APPLIED FIELD 

Suppose the applied magnetic field has a small trans- 
verse component, so that Ha = Hix + H2y, where 



H2 — eh2{x) 



(52) 



(/i2 depends on x but not t). For simplicity, let K2 = 
0. Let m.e{x,t) denote the solution of the LL equation 
with initial condition mg(x, 0) — m*(x). As above, let 
T{xe{t))R{<j)e{t))T[n^ denote the translated and rotated 
optimal profile closest to nie at time t. 

We first note that, unless /12 vanishes as x ±00, 
nie will not remain close to T{xe{t))R{(j)e{t))m^,. For 
example, if /12 is constant, then away from the domain 
wall, nig will relax to one of the local minimizers of the 
homogeneous energy Ki{l — mf) — Ha • m, and these 
do not lie along ±x for H2 7^ 0. It follows that ||m£ — 
T{xe{t))R{4>i{t))In.^,\\ will diverge with time. 

Physically, this divergence is spurious. It stems from 
the fact that we are taking the wire to be of infinite ex- 
tent. One way to resolve the issue, of course, would be to 
take the wire to be of finite length. However, one would 
then no longer have an explicit analytic solution of the 
LL equation. 



Here we shall take a simpler approach, and assume 
that the transverse field h2{x) approaches zero as x ap- 
proaches ±00. In fact, for technical reasons, it will be 
convenient to assume that the integral of /i2-l-/i2^, i.e. the 
squared Sobolev norm ||/i2||ffi, is finite. Then without 
loss of generality, we may assume 



(hl + h'2')d^ = l. 



(53) 



Under this assumption, the main result of this section is 
that nie stays close to an optimal profile up to translation 
and rotation. That is, for some Ci > 0, 



(54) 



The demonstration proceeds as in the preceding sec- 
tion, so we will discuss only the points at which the 
present case is different. The main difference is that, 
in place of Eq. we get (by considering the LL equa- 
tion with H2 ^ rather than K2 ^ 0) the following 
expression for G{m^) to leading order O(e^): 

G{m,) = e'^[acos(t)4t) J cosO^ {6'( - cos26.Ji)h2dS, 

— Q;sin0*(<) j sin0* {(!)'{ — 2 cos9^(j)'i) h2 dS, 

-sin^^t) J {e'l - cos2eji)h2d^ 

— cos(/)*(t) J sin 6^ COS0* (0" — 2 cos ^^t/i'^) /i2 . 

(55) 

After some straightforward manipulations including in- 
tegration by parts, and making use of the inequality 
Eq. ([50)1 , one can show that 



cos 0* 



cos 26 ^9i) /i2 c?^ 



< 



/3, 



1 

2/3' 



sin 9^, {(j>i — 2 cos 6** (p'l ) /12 

/3„ 



. 

< - sm ( 
- 2' 



1 

2/3' 



cos 20,6*1) h2 d^ 



< 



1 

2^' 



sin 9^ cos 9^ (0" — 2 cos 0, ) /12 



< ^\\sin9,^[\\^ 



1 

2/3 



TTa- (56) 



From Eqs. ^ and ^ it follows that 



and 



+ sin^ 9,(j)'^^^ d^<^AE„ 



(57) 



(58) 
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Substituting Eqs. (l56])-(l58l) into Eq. dSS]), we get that 

|G(m,)| < (1 + a) (^3/3Ai?, + ie2^ . (59) 

This estimate is of the same form as and the ar- 

gument given there, with (3 chosen appropriately, estab- 
hshes Eq. (|54l) . 

V. NUMERICAL STUDIES 

In the preceding Sections HIHlVI we have shown that 
the precessing solution is linearly stable; to leading or- 
der 0(e), a perturbed solution either approaches or stays 
close to the precessing solution up to a translation and 
rotation, according to whether the perturbation is to the 
initial conditions or to the anistropy and transverse ap- 
plied magnetic field in the LL equation. Here we present 
numerical results which verify nonlinear stability for the 
precessing solution under small perturbations. To this 
end, we investigate the energy, AE^ — E{m^)~E{ni^), of 
the numerically computed perturbed DW, m^{x,t), rela- 
tive to the minimum energy i?(m*) of an optimal profile, 
as a function of time t. Throughout, E is taken to be 
the uniaxial micromagnetic energy given by Eq. ([1} with 
K2 = 0. As in the preceding sections, we choose units so 
that A = Ki = 1. In these units, £^(m*) = 2. In typical 
ferromagnetic microstructures, the value of the Gilbert 
damping parameter a is known to lie between 0.04 and 
0.22 (see e.g. RefJ^ and references within), so we take 
a = 0.1 throughout our numerical study. 

A. Perturbed initial profile 

We first investigate the evolution of a DW, m^{x,t), 
from an initial perturbation of an optimal profile. We 
take the initial condition in polar coordinates to be given 
by 

e,{x,0) = (^-^-^^ , Mx)^c^o + e2X, (60) 

which corresponds to stretching the unperturbed profile 
along and twisting it around the axis of the nanowire. 
The applied field is directed along the nanowire. Ha = 
ffix, and we take K2 = 0. 

Figure [1] shows the dependence of the relative energy 
AEf: on time t for different values of the applied field Hi . 
The figure presents 13 curves corresponding, from top to 
bottom, to Hi varying from —1.2 to at the increment 
of 0.1. In the initial condition given by Eq. (|60| . we take 
ei = 0.1 and €2 = 7r/50. 

Figure [1] clearly indicates that AEf{t) decays exponen- 
tially for weak applied fields, \Hi \ < 1/2, in accord with 
the analytic result Eq. ([311 . However, for \Hi \ ~ 1, devi- 
ations from exponential decay are evident, and the pre- 
cessing solution appears to become unstable for \Hi | ^ 1. 




50 100 150 200 

t 



FIG. 1: (Color online) Relative energy, AE^{t), of the per- 
turbed DW for 13 different values of the applied field Hi . See 
text for discussion. 

B. Small hard-cixis anisotropy 

We consider next the evolution of a DW from an opti- 
mal profile at i = when the hard-axis anisotropy K2 is 
nonvanishing. We fix Hi = —0.5. 




FIG. 2; (Color online) Relative energy, AEt{t), of the per- 
turbed DW for 5 different values of the hard-axis anisotropy 
constant K2- See text for discussion. 

Figure [2] shows the dependence of the relative energy 
AEf: on time t for different values of K2. The figure 
presents 5 curves corresponding, from top to bottom, to 
K2 varying from 0.1 to 0.02 at the decrement of 0.02. 
(The blue and red colorings alternate to make adjacent 
curves more easily distinguishable.) It is evident that the 
relative energy remains small, verifying the linear analy- 
sis of Sec. [ml 

Figure [31 shows the maximum value of the relative en- 
ergy AEf (over the interval < t < 80) as a func- 
tion of K2 . Red squares represent numerically computed 
values. The black solid curve is the parabola Cif/f|, 
with Ck — 1.3207 fitted by the method of least squares 
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FIG. 3: (Color online) Maximum value of the relative en- 
ergy Ai5e of the perturbed DW as a function of the hard- 
axis anisotropy K2- Numerically computed values are repre- 
sented by (red) squares. The (black) solid curve is a parabola, 
max(A£;s) = CkKI with Ck = 1.3207, fitted by the method 
of least squares through the data points with K2 < 0.04. 

through the data points with K2 < 0.04. We obtain con- 
vincing confirmation of the leading-order analytical re- 
sult Eq. For larger values of K2, we see departures 
from quadratic dependence; for sufficiently large values 
of if 2 (not shown), the Walker solution was recovered. 



C. Small transverse applied field 

Finally, we address the stability of the precessing solu- 
tion under an applied magnetic field. Ha — Hix + H2y, 
with a small transverse component, H2{x). As discussed 
in Sec. IIVI we want H2{x) to vanish as a; — >■ ±00. Here 
we take 

H2{x) ^ H2w{x), (61) 

where w{x) is equal to one inside the window < a; < 20 
and vanishes outside (the argument of Section lTVl is easily 
modified to establish the linear stability result Eq. p8)) 
in this case). We consider the evolution of a DW given 
at t = by the optimal profile m, centred at a; = 0. We 
take Hi = —0.5, so that in the absence of the transverse 
field, the DW velocity is positive (cf. Eq. (|13l) ) and the 
DW crosses the window. We take K2 — 0. 

Figure |3] shows the dependence of the relative energy 
AE'g on time t for different values of the transverse field 
amplitude H2. The figure presents 5 curves correspond- 
ing, from top to bottom, to H2 varying from 0.1 to 0.02 
at the decrement of 0.02. (The blue and red colorings al- 
ternate to make adjacent curves more easily distinguish- 
able.) The relative energy AE^{t) is presented over the 
time interval < t < 400, which, for small values of 
H2, is sufficient for the DW to traverse the spatial win- 
dow < X < 20 (cf. Eq. (HH)). The results confirm that 




t 



FIG. 4; (Color online) Relative energy, AE^{t), of the per- 
turbed DW for 5 different values of the transverse field am- 
plitude H2- See text for discussion. 



the relative energy of the perturbed magnetization pro- 
file remains small for small values of H2, in accord with 
the leading-order results of Section HVl 




S2 



FIG. 5: (Color online) Maximum value of the relative energy 
AEe of the perturbed DW as a function of the amplitude 
of the transverse applied field, H2- Numerically computed 
values are represented by (red) squares. The (black) solid 
curve is a parabola, meix{AE^) = ChHI with Ch = 99.6586, 
fitted by the method of least squares through the data points 
with H2 < 0.04. 

Figure [S] shows the maximum value of the relative en- 
ergy Ai?e (over the interval < t < 400) as a function of 
H2. Red squares represent numerically computed values. 
The black solid curve corresponds to the parabola C////| 
with Ch — 99.6586 fitted by the method of least squares 
through the data points with H2 < 0.04. The figure 
provides a confirmation of the leading-order analytical 
result of Sec. IIVI that the maximum relative energy de- 
pends quadratically on H2 for small H2 ■ Deviations from 
the parabolic dependence can be seen for H2 > 0.08. 
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VI. CONCLUSIONS 

The precessing solution is a new, recently reported ex- 
act solution of the Landau-Lifschitz-Gilbert equation. It 
describes the evolution of a magnetic domain wall in a 
one-dimensional wire with uniaxial anisotropy subject to 
a spatially uniform but time-varying applied magnetic 
field along the wire. We have analysed the stability of the 
precessing solution. We have proved linear stability with 
respect to small perturbations of the initial conditions as 



well as to small hard- axis anisotropy and small transverse 
applied fields, provided the applied magnetic field along 
the wire is not too large. We have also carried out nu- 
merical calculations that confirm full nonlinear stability 
under these perturbations. 

Numerical calculations suggest that, for sufficiently 

large perturbations and applied longitudinal fields, the 
precessing solution becomes unstable, and new stable so- 
lutions appear. It would be interesting to analyse these 
bifurcations and study these new regimes for DW motion. 
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